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Abstract. It is shown that the traditionally accepted "Rouse values" for the critical 
exponents at the gelation transition do not arise from the Rouse model for gelling 
polymers. The true critical behaviour of the Rouse model for gelling polymers is 
obtained from spectral properties of the connectivity matrix of the fractal clusters 
that are formed by the molecules. The required spectral properties are related to the 
return probability of a "blind ant" -random walk on the critical percolating cluster. The 
resulting scaling relations express the critical exponents of the shear-stress-relaxation 
function, and hence those of the shear viscosity and of the first normal stress coefficient, 
in terms of the spectral dimension d s of the critical percolating cluster and the 
exponents a and r of the cluster-size distribution. 

PACS numbers: 64.60.Ak, 64.60.Ht 
1. Introduction and result 

During a chemical gelation reaction the crosslinking of (macro-) molecules leads to the 
formation of permanent clusters of very different shapes and sizes. In the absence of a 
macroscopically large cluster such a system acts like a liquid (sol). Otherwise one deals 
with an amorphous solid state (gel). Hence, the sol-gel or gelation transition appears 
as a percolation transition [1-4], which is driven by crosslink concentration. 

Gelling liquids exhibit striking mechanical properties which have been continuously 
studied over the years by experiments [5-12], theories [13-24] and simulations [25,26]. 
For example, when subjected to a homogeneous shear flow, distinct relaxation patterns 
are observed, which are due to the participation of many different excitation modes. 
More precisely, experiments suggest the scaling form [2,4,6-8, 12, 16] 

G{t) ~r A g(t/t*) with t*{e)~ £ - z (1) 

for the (shear-) stress-relaxation function in the sol phase for asymptotically long times t 
and crosslink concentrations c close to the critical point c cr j t , i.e. for e := (c cr i t — c)/c cr j t <C 
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1. The typical relaxation time t* diverges with a critical exponent z > for e J, 0. The 
scaling function g is of order unity for small arguments so that one finds the algebraic 
decay G(t) ~ t~ A with a critical exponent A > for t — > oo at the critical point 
£ = 0. For large arguments g decreases faster than any inverse power. Sometimes a 
stretched exponential has been proposed [8, 16] for g, implying the long-time asymptotics 
- lnG(t) ~ (t/t*) 7 with some < 7 < 1 in the sol phase. 

Within linear viscoelasticity one finds 77 = f™dtG(t) for the zero-frequency 
viscosity and ty± = 2 f™dttG(t) for the first normal stress coefficient [4, 27]. The scaling 
form (1) of the stress- relaxation function determines their critical divergencies rj ~ s~ k 
and ^1 ~ e~ e as e [ by the scaling relations [4, 24] 

k — z(l — A) and £ = z(2 - A) = k + z . (2) 

Thus, it suffices to know any two of the four critical exponents A, z, k and I. 

The Rouse model is one of the most basic models for polymer dynamics. By 
definition [28-30], it is a bead-spring model for linear polymer chains with no other 
interactions among monomers (i.e. beads) than connectivity along the chain, as 
accounted for by harmonic springs. Moreover, each bead behaves as a Brownian particle 
in that it experiences a friction force and a thermal noise force. Despite its simplicity 
the Rouse model has provided a much valued description of linear polymer chains over 
the last fifty years. Its generalization from linear chains to arbitrarily connected clusters 
of monomers is straightforward and will be recalled in Section 2. Static and dynamic 
properties of such Rouse clusters have been discussed, for example, in [15,31-34], and 
recent applications to gelling polymer melts can be found in [18,20-24]. More than 
twenty years earlier, de Gennes was the first who aimed at a description of the gelation 
transition in terms of the Rouse model. Guided by scaling arguments, he claimed the 
relation k = 2v — f3 for the critical exponent of the viscosity [14]. Here, the exponent 
v governs the divergence of the correlation length and (3 the growth of the gel fraction. 
Using this result, it was also argued [6, 7, 10, 16, 35] that A = dvj (dv + k), where d is the 
space dimension. Inserting the numerical values for three-dimensional bond percolation 
[36,37], as is generally accepted, gives k pa 1.35 and A pa 0.66. These results have 
become known as the "Rouse values" for the critical exponents and are commonly 
quoted in the literature, see e.g. the reviews [1-4]. 

In contrast, the exact solution of the Rouse model for gelling polymers in [20,21] 
leads to 

k = (l-r + 2/d s )/a (3) 

and thus disproves the earlier scaling arguments — and hence the Rouse value — for k 
mentioned above. Here d s denotes the spectral dimension [36-38] of the incipient 
percolating cluster, and the exponents a and r are related to the cluster-size distribution, 
as defined in (12) below. The derivation of (3) in [20,21] exploits a correspondence 
between Rouse clusters and an electrical resistor network, which rests on the linearity 
of Hooke's law and Ohm's law. This correspondence is exact. It differs from de 
Gennes' superconductor/conductor- network analogy [39], as can be seen most clearly 
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Table 1. Critical exponents for a gelling Rouse melt based on (4) and (2). 



A = d s (r-l)/2 z = 2/(d s a) k = z(l — A) £ = z(2 - A) 



2D 
3D 



0.70 
0.79 
1 



3.8 
3.3 
3 



1.2 

0.71 





5.0 
4.1 

3 



MF 



by comparing the resulting numerical values for k in two space dimensions [20]. 
Unfortunately, the approach of [20,21], which led to (3), does not seem to be suitable 
for a computation of the other exponents A, z or £ — except for the case where the 
statistical ensemble of the clusters corresponds to mean-field percolation [22-24]. Hence, 
the authors of [24] had to use the numerically estimated value A m 0.83, which was 
obtained in [22] for the Rouse model for gelling polymers with clusters distributed 
according to three-dimensional percolation. 

The present Paper fills this gap. We present a simple analytical calculation which 
determines both A and z — and hence also k and £ by (2) — for Rouse clusters distributed 
according to <i-dimensional percolation. This will be achieved by a well-known relation 
between the dynamics of Rouse clusters and a "blind ant" -random walk on percolation 
clusters [15,36,37]. In particular, we will express the stress-relaxation function G in 
terms of the return probability of such a random walk. In this way we will deduce the 
scaling relations 



for the microscopic Rouse model for gelling polymers. Of course, (4) and (2) are 
consistent with existing exact results, that is, equation (3) for k [20, 21] and the mean- 
field values for A, z and £ [22-24]. Table 1 summarizes the exact values of the critical 
exponents pertaining to the Rouse model for gelling polymers with clusters distributed 
according to two- and three-dimensional lattice-bond percolation, respectively mean- 
field percolation. 

The result for A in (4) differs from the relation A = du/(du + k) that has grown 
out of the scaling arguments mentioned above. Hence, it is not only the "Rouse value" 
for k but also the one for A that has been ascribed incorrectly to the Rouse model for 
gelling polymers. 

The rest of the Paper is organized as follows. After recalling the Rouse model for 
gelling polymers in Section 2, we will derive the scaling relations (4) within this model 
in Sections 3 to 6. Section 7 contains a discussion and an outlook. 

2. Rouse model for gelling polymers 




and 



2 
d s a 



(4) 



In order to set up the notation we recall the standard definition [15, 29, 30, 32-34] of the 
Rouse model for clusters of monomers with an arbitrary topology. Previous applications 
to gelling polymer melts can be found in [18, 20-24]. We consider N point-like monomers 
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o = -C 



^Ri(t)-v ext (Ri(t),t) 



<t)+Ut) (6) 



with positions Ri(t), i = 1,...,N, in d-dimensional space, d > 2. Their motion is 
constrained by M randomly chosen, permanent, harmonic crosslinks which connect the 
pairs of monomers (i e ,j e ), e = 1, . . . , M, and give rise to the potential energy 

, M , N 

U '-=^ " Rjf =: ^2 E r « * • * ■ ( 5 ) 

e=l i,j 

Here, we have chosen units such that the inverse temperature (k^T)" 1 is equal to one, 
and the fixed length a > serves as an inverse crosslink strength. The second equality 
in (5) introduces the random N x iV-connectivity matrix T, which encodes all properties 
of a given crosslink realisation. The simple relaxational dynamics 

dU 
dRi 

governs the motion of the monomers in the externally applied shear flow v^ xt (r,t) : = 
K>(t) S a ,i r 2, which is oriented along the 1-direction and involves the spatially 
homogeneous, time-dependent shear rate K,(t). Here, Greek indices label Cartesian 
components, r = (ri,...,^) and 5 a ^ denotes the Kronecker symbol. A friction 
force with friction constant ( applies to monomer i, if its velocity deviates from 
that of the surrounding shear flow. Moreover, the crosslinks exert the Hookean 
forces —dU/dRi, and the last contribution in (6) is due to a randomly fluctuating 
thermal-noise force £i(t), obeying Gaussian statistics with zero mean and covariance 
£°(t)£j(t') = 2(5 a ^5ij5(t — t'). Here, 5(t) denotes the Dirac delta function. Given the 
homogeneous shear flow v cxt and any fixed configuration C := {(i e , j e )}^L l of permanent 
crosslinks, the equation of motion (6) is linear and can be solved exactly for every 
realisation of the thermal noise [21]. 

It remains to specify the statistical ensemble which determines the probabilistic 
occurrence of different crosslink configurations C. Here we consider any bond- 
percolation ensemble which is amenable to a scaling description, does not distinguish 
any site (i.e. monomer), yields a fractal percolating cluster at criticality and where the 
maximum number of bonds (i.e. crosslinks) emanating from a given site is uniformly 
bounded. Besides certain versions of continuum percolation [40], this also includes 
bond percolation on the simple cubic lattice in d dimensions and on the Cayley tree 
[36,37]. The statistical average over all crosslink configurations will be denoted by 
angular brackets ( • ) and, when using this notation, we implicitly assume that the 
macroscopic limit iV — > oo, M — > oo, M/N — > c is carried out, too. 



3. Stress relaxation in the Rouse model for gelling polymers 

The crosslinks exert shear stress on the molecular clusters due to the externally applied 
shear flow. Following Kirkwood [29,30], the Cartesian components of the stress tensor 
are given in terms of a force-position correlation 
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Here, po denotes the density of monomers, and one has to insert the explicitly known 
[21] solutions Ri(t) of the Rouse equation (6) with some initial data Ri(t ). Finally, 
the limit to ^ —oo ensures that all transient effects stemming from the initial data will 
disappear. This yields the exact result 

( 2 fidsK(s) 1 ••• \ 

(a(t)} = G(0)1+ / dt' G(t - t') /c(t') 

J — oo 



1 











V 



(8) 



/ 



for the crosslink-averaged, macroscopic stress tensor, where 1 is the 3 x 3-unit matrix 
and 

6t 



G(t) = G(t,e) 



Po 
N 



Tr 



1-K 



o) exp 



Co? 



(9) 



represents the stress-relaxation function [21,24]. Here, Tr is the trace and E denotes 
the eigenprojector on the space of zero eigenvalues of T, which correspond to translations 
of whole clusters [18, 20, 21]. Within the simple Rouse model the zero eigenvalues do not 
contribute to stress relaxation because there is no force acting between different clusters. 
The only contribution to stress relaxation is due to deformations of the clusters. 



4. Cluster decomposition 

A crosslink configuration C can be decomposed into K = K(C) disjoint clusters, and 
we denote the number of monomers in its fcth cluster Af k = A4(C) by N k = N k {C). 
The connectivity matrix V(C) of the whole configuration is thus a block- diagonal sum of 
connectivity matrices T(Af k ) of different clusters. Hence, when evaluating the trace 
in (9), one gets an additive contribution from every cluster. Upon sorting these 
contributions with respect to the number of monomers in the clusters and introducing 
the mean number density r n := (iV _1 Y2k=i ^N k ,n) of clusters with n monomers — also 
known as cluster-size distribution — , we arrive at the cluster-decomposed form 

oo 

G(t,e) = J2 n ^(e)G n (t,e) (10) 

n=l 



G n (t,e) :=(^Tr 



E Q ) exp ( -— r 



(11) 



of the stress-relaxation function. Here, the mean stress-relaxation function of n-clusters 

6t 

is defined in terms of the (conditional) average (A) n := t^ 1 (N~ 1 Y^k=i $N k ,nA(J\fk)) for 
an observable A over all clusters with n monomers. The cluster-size distribution r n gives 
rise to the critical exponents a and r through its scaling form 

T n (e)~n- T f\(n/n*(e)) (12) 

for n 3> 1 and crosslink concentrations close to the critical point [36,37]. It involves 
the typical cluster mass n*(e) ~ e~ l l a ', which diverges for £ I 0. The scaling function fi 
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decreases faster than any inverse power for large arguments, whereas for small arguments 
it approaches a non-zero, finite constant. It will turn out later that in order to calculate 
the critical exponents A and z of G, it suffices to know G n at the gelation transition. 
For this reason we will only be concerned with critical percolation clusters in the next 
paragraph. 



5. Blind ant on critical percolation clusters 

Consider a random walker — coined "blind ant" by de Gennes [41] — that moves along 
a bond from one site to another in the same cluster at discrete time steps [36-38,42]. 
More precisely, let b denote the maximum number of possible bonds which are allowed 
to emanate from a site in the percolation model under consideration. If the ant happens 
to visit site % at time s, which is connected with bi < b bonds to other sites, then it will 
move with equal probability 1/6 along any one of the bi bonds within the next time step 
and stay at site i with probability 1 — bi/b. By the definition (5) of the connectivity 
matrix T of the cluster, one has Tu = bi for its diagonal matrix elements, T^- = — 1 if two 
different sites i ^ j are connected by a bond and zero otherwise. Hence, the associated 
master equation for the ant's sojourn probability Pi(s) for site i at time s reads 

Pt (s + 1) = (1 - r it /b) Pl (s) + Y^i-Vij/bMs) , (13) 

which is equivalent to 

Pl (s + 1) - Pt (s) = -b- 1 J2 FijPjis) . (14) 

3 

Here the summation extends over all sites in the cluster. For long times s 1, it 
is legitimate to replace the difference (quotient) on the left-hand side of (14) by a 
derivative. This yields the solution Pi(s) = le~ sT ^ b ]u , which corresponds to the initial 
condition Pi(Q) = Si jio . 

Next we consider P ( - n \s) := (pi (s)) n \ £=0 , the mean return probability to the 
starting point after time s, where the average is taken over all critical percolation clusters 
with n sites. Likewise, p(°°\s) stands for the mean return probability on the incipient 
percolating cluster. Clearly, these definitions are independent of the starting point io, 
because on average there is no distinguished site by assumption. Thus we can also write 

P (n) (s) = /-Tre" sr/ A (15) 

\ U I n e =0 

for finite n. The decay of P(°°\s) for long times s 3> 1 is determined by the inverse 
number of distinct sites which the ant has visited up to time s and, thus, by the spectral 
dimension d s of the incipient percolating cluster according to [37,38,42] 

P(°°)(s) ~s- d °/ 2 . (16) 

The fractal nature of percolation clusters at the critical point implies that clusters 
with n»l sites typically look like regions of n connected sites that are cut out of the 
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infinite cluster. Therefore we conclude P( n \s) ~ s d ^ 2 for times 1 < s < s„, where 
[37, 38,42] 

s n ~ n 2/d ° (17) 

is the crossover time at which the ant first recognizes that it does not move on the 
infinite cluster, but on a finite one with n sites, i.e. the time it needs to travel a distance 
that is comparable to the cluster's extension. For times s ^> s n , on the other hand, 
every site in the cluster is equally likely to be visited by the ant, hence p( n \s) ~ 1/n 
in this case. Taken together, we can write 

P (n \s)~s- d °/ 2 f 2 (s/s n ) + l/n (18) 

for s >> 1. The cut-off function /b^c) is of order one for x < 1 and decreases rapidly to 
zero for x — * oo. Upon comparing (15) and (18), we infer the relation 

~ s- ds/2 f2(s/s n ) (19) 

£ = 

for s 1. To derive (19) we have also used the explicit form (-E'o)y — n' 1 , 1 < i,j < n, 
of the matrix elements of the eigenprojector corresponding to the non-degenerate zero 
eigenvalue of V for the case of a single cluster with n sites. 



-Tr[(l-£ )e- sr / b ] 
n 



6. Critical behaviour of stress relaxation 



The long-time decay of the stress relaxation function G is determined by the biggest 
clusters. Close to the gelation transition we may thus use the scaling form (12) of the 
cluster-size distribution and approximate the (Riemann) sum on the right-hand side of 
(10) by the corresponding integral 

poo 

G(t,e)~(n*(e)) 2 ^ dxx^f^G^^e). (20) 
Jo 

If (20) is to be consistent with the scaling form (1), we are forced to require the functional 
dependence 

G n (t,e) ~t- a f 3 (n/n*(e),t/t*(E)). (21) 

The yet unknown exponent a and the time scale t*(e) ~ e~ z can be determined from 
the limiting case e = 0, 

G n (t, 0) ~ t- d *l 2 f 2 (t/t n ) with t n ~ n 2 ' d ° , (22) 

which is obtained from (11), (19) and (17). Indeed, for (22) to be the limit of (21) as 
e I 0, one has to identify a = d s /2, and the diverging e-dependencies in / 3 have to cancel 
each other in order to yield a non-trivial limit. This is only possible if fs(x, y) ~ $2(31/ x a ) 
as x, y — > with some exponent a which guarantees t*(e) ~ (n*(e)) a . Consequently, we 
also get the relation t n ~ n a . Hence, a = 2/d s by (22), and we conclude t*(e) ~ s~ z 
with z = 2/(d s a). Moreover, expressing n*(e) in terms of t*(e), the stress- relaxation 
function (20) takes on the form 

POO 

Git, e) ~ r A (t/t*(e)) A - d °/ 2 / dx x^hix) f 3 (x, t/t*(e)) (23) 

Jo 
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for e <C 1 and t — > oo with A = (t — l)d s /2. This completes the derivation of (4). 
7. Discussion and outlook 

Table 1 summarizes the exact critical exponents for stress relaxation in the Rouse model 
for gelling polymers below and at the gel point. The exponents differ from the predictions 
k = 2v — (3 and A = dvj [dv + k) of earlier scaling arguments [6, 7, 10, 14, 16, 35]. What 
is the reason for this discrepancy? The scaling arguments involve the Hausdorff fractal 
dimension df := d—(3/v of rigid percolation clusters at c cr i t . Rouse clusters, however, are 
thermally stabilized, Gaussian phantom clusters with the Hausdorff fractal dimension 
df := 2d s /{2 — d s ) [15,31,33], which is different from df in space dimensions below 

(G) 

the upper critical dimension d u = 6. Indeed, if one replaces df by df in these scaling 
arguments, as one should consistently do in a Rouse description, the results will coincide 
with Table 1. 

Since the scaling relations k = 2v — (3 and A = dvj (du + k) involve the Hausdorff 
fractal dimension df of rigid percolation clusters, it is sometimes argued that they 
describe the behaviour of a more realistic model, which, in addition to the interactions 
of the Rouse model, accounts for excluded- volume effects, too, see e.g. [35]. To the 
best of the author's knowledge this claim has not been verified yet by strict analytical 
arguments within a microscopic model. One may even have doubts whether this claim is 
generally true: Extensive molecular-dynamics simulations [26] of a system of crosslinked 
soft spheres in three dimensions, with a cluster ensemble as in the present Paper and an 
additional strongly repulsive interaction at short distances, yield the values k ~ 0.7 and 
A 0.75, which are remarkably close to the predictions of the Rouse model for gelling 
polymers, see Table 1. On the other hand, simulations of the bond-fluctuation model 
in [25] imply k ~ 1.3 and are thus in favour of the claim. However, the viscosity is 
not measured directly in these latter simulations. Rather it is derived from the scaling 
of diffusion constants and an additional scaling assumption that may be questioned 
[26] . Hence, it is an open problem to what extent the critical exponents of Table 1 are 
modified by excluded-volume interactions. 
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